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HIGHLIGHTS 


►  A  dynamic  electrochemical  polarization  battery  model  is  used  for  state  estimation. 

►  An  improved  parameter  identification  test  is  used  to  get  accurate  model’s  parameters. 

►  The  continuous  peak  power  capability  estimation  approach  is  proposed. 

►  An  AEKF-based  SoC  and  peak  power  capability  joint  estimation  approach  is  proposed. 

►  The  robustness  and  reliability  of  the  six  steps  joint  estimator  are  evaluated. 
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This  paper  uses  an  adaptive  extended  Kalman  filter  (AEKF)-based  method  to  jointly  estimate  the  State  of 
Charge  (SoC)  and  peak  power  capability  of  a  lithium-ion  battery  in  plug-in  hybrid  electric  vehicles 
(PHEVs).  First,  to  strengthen  the  links  of  the  model’s  performance  with  battery's  SoC,  a  dynamic  elec¬ 
trochemical  polarization  battery  model  is  employed  for  the  state  estimations.  To  get  accurate  parameters, 
we  use  four  different  charge-discharge  current  to  improve  the  hybrid  power  pulse  characteristic  test. 
Second,  the  AEKF-based  method  is  employed  to  achieve  a  robust  SoC  estimation.  Third,  due  to  the 
PHEVs  require  continuous  peak  power  for  acceleration,  regenerative  braking  and  gradient  climbing,  the 
continuous  peak  power  capability  estimation  approach  is  proposed.  And  to  improve  its  applicability, 
a  general  framework  for  six-step  joint  estimation  approach  for  SoC  and  peak  power  capability  is 
proposed.  Lastly,  a  dynamic  cycle  test  based  on  the  urban  dynamometer  driving  schedule  is  performed  to 
evaluate  the  real-time  performance  and  robustness  of  the  joint  estimation  approach.  The  results  show 
that  the  proposed  approach  can  not  only  achieve  an  accurate  SoC  estimate  and  its  estimation  error  is 
below  0.02  especially  with  big  initial  SoC  error;  but  also  gives  reliable  and  robust  peak  power  capability 
estimate. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Energy  crisis,  environmental  issues  have  promoted  the  devel¬ 
opment  of  various  types  of  new  energy  vehicles,  which  has  been 
listed  as  one  of  the  seven  strategic  emerging  industries  in  China 
[1,2],  Plug-in  hybrid  electric  vehicles  (PHEVs)  are  emerging  as 
important  personal  transportation  options  for  petroleum  displace¬ 
ment  and  energy  diversification,  and  its  double  operation  modes, 
hybrid  electric  vehicle  mode  for  suburbs  or  other  driving  cycles  and 
pure  electric  vehicle  mode  for  urban,  is  favorable;  to  meet  these 
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requirements,  the  battery  system  in  PHEVs  requires  both  high 
specific  power  and  high  energy  [3-5]. 

The  lithium-ion  battery  is  widely  used  in  many  fields  because  of 
its  advantages  of  high  voltage,  long  cycle-life,  high  specific  energy 
and  high  specific  power,  which  also  makes  it  a  promising  candidate 
for  the  PHEVs  [6,7],  For  PHEVs,  on  the  one  hand,  battery  State  of 
Charge  (SoC)  is  an  important  quantity,  as  it  is  a  gauge  of  the 
remaining  capacity  in  the  battery.  Thus,  in  order  to  manage  the 
battery  efficiently,  an  accurate  SoC  estimation  method  is  one 
of  the  key  issues.  In  addition,  the  accurate  SoC  estimation  can 
improve  the  power  distribution  efficiency  greatly,  and  also  is  the 
premise  of  the  accurate  estimates  of  the  battery  peak  power 
capability  [8—12],  On  the  other  hand,  accurate  peak  power  esti¬ 
mates  are  critical  in  practical  applications  since  it  is  necessary  to 
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determine  the  power  available  in  the  moment  to  meet  the  accel¬ 
eration,  regenerative  braking  and  gradient  climbing  power 
requirements  without  fear  of  over-charging  or  over-discharging  the 
battery.  More  importantly,  accurate  peak  power  capability  esti¬ 
mates  for  the  battery  will  optimize  the  relation  between  the  battery 
capacity  and  size  with  vehicle’s  performance,  which  benefits  the 
vehicle’s  general  potency  [2,13,14],  Therefore,  to  provide  an  effi¬ 
cient  guarantee  for  the  optimization  of  PHEVs  power  systems,  an 
accurate  online  SoC  and  peak  power  capability  calculation  algo¬ 
rithm  is  particularly  critical  to  maintain  optimum  battery 
performance. 

However,  due  to  the  strong  time-variable  and  nonlinear  perfor¬ 
mance  of  battery  system,  it  is  hard  to  measure  its  SoC  directly  for  the 
various  influence  factors  and  complicated  electrochemical  process 
from  the  practice  application.  Some  SoC  estimation  approaches 
which  are  presented  with  the  development  of  battery  powered 
electric  vehicles  can  be  generally  classified  into  three  kinds.  The  first 
is  the  direct  measurement  method,  which  involves  measurements  of 
voltage  and  electrochemical  impedance  [15,16],  However,  due  to  the 
uncertain  driving  cycles,  complex  application  environments,  it  is 
hard  to  measure  the  characterization  parameters  real-timely,  and 
there  is  much  to  be  solved  for  practical  application.  The  second  is 
ampere-hour  counting  method,  which  is  based  on  both  current 
measurement  and  integration  [17],  However,  its  performance  is 
highly  dependent  on  the  measuring  accuracy  of  current,  and  this 
open-loop  calculation  method  can  easily  lead  to  accumulated 
calculation  errors  due  to  uncertain  disturbances  from  the  practice 
application  and  lack  of  necessary  corrective  resolution.  The  third  is 
the  model-based  method  with  filter  algorithms,  which  carry  out 
estimation  by  means  of  state-space  battery  models  [18—24],  In  Refs. 
[8—10],  the  authors  use  the  extended  Kalman  filter  (EKF)  to  adap¬ 
tively  estimate  the  SoC  based  on  a  simplified  battery  models. 
However,  the  Kalman  filters-based  algorithm  strongly  depends  on 
the  precision  of  the  battery  model  and  the  predetermined  variables 
of  the  system  noise  such  as  mean  value,  relevance  and  covariance 
matrix.  An  inappropriate  information  matrix  of  the  system  noise 
may  lead  to  remarkable  errors  and  divergence  [12],  Therefore,  an 
adaptive  extended  Kalman  filter  (AEKF)-based  method  has  been 
applied  to  implement  online  SoC  estimation  in  Refs.  [12,23]  to 
improve  the  estimate  accuracy  by  adaptively  updating  the  process 
and  measurement  noise  covariance. 

On  the  other  hand,  some  peak  power  capability  estimation 
approaches  are  presented  with  the  development  of  electric  vehicles 
technology  [2,13,14],  The  most  commonly  used  approach  is  hybrid 
pulse  power  characterization  (HPPC)  method  proposed  by  the  Idaho 
National  Engineering  &  Environmental  Laboratory,  which  deter¬ 
mines  the  static  peak  power  in  laboratory  environments,  and  the 
estimates  are  over  optimism.  However,  it  is  not  suitable  to  estimate 
the  continuous  peak  currents  that  are  available  for  the  next 
sampling  interval  At,  additionally,  the  method  neglects  design 
limits  like  cell  current,  cell  power  or  SoC  [2],  As  an  improvement, 
the  voltage-limited  method  was  proposed  by  Plett  [14],  However, 
these  two  Rint  model-based  methods  can  hardly  simulate  the 
dynamic  voltage  performance  and  the  estimates  will  diverge  from 
the  practical  capability.  To  solve  this  problem,  the  authors  [13] 
propose  a  dynamic  electrochemical  polarization  (EP)  battery 
model-based  multi-parameter  peak  power  capability  estimation 
method,  and  the  experiment  and  simulation  results  show  that  the 
method  has  better  performance  than  other  commonly  used 
methods.  However,  this  method  can  efficiently  estimate  the  peak 
power  capability  of  the  battery  for  the  next  sampling  interval  At, 
the  estimates  of  multiple  sampling  intervals  Ats  can  hardly  be 
achieved,  which  is  because  the  battery  terminal  voltage  will  change 
at  any  sampling  interval  when  battery  is  discharged  or  charged  at 
any  current  or  power,  therefore  it  is  an  inappropriate  way  to  replace 


the  sampling  interval  At  with  multiple  sampling  intervals  Ats  to  get 
a  longer  time  estimation.  As  a  result,  this  method  fails  to  provide 
continuous  power  for  multiple  sampling  intervals,  such  as  15  s,  30  s, 
60  s  etc.  when  acceleration,  regenerative  braking  or  gradient 
climbing. 

In  this  paper,  based  on  the  dynamic  electrochemical  polariza¬ 
tion  (EP)  battery  model,  an  AEKF-based  SoC  and  peak  power 
capability  joint  estimation  approach  for  a  3.7  V/35  Ah  LiMn204 
lithium-ion  battery  used  in  PHEVs  is  proposed.  A  description  of  the 
EP  battery  model  and  its  parameters  identification  process  are 
given  in  Section  2.  The  six  steps  EP  model-based  SoC  and  peak 
power  capability  joint  estimation  approach  with  AEKF  algorithm  is 
depicted  in  Section  3.  The  experiment  and  evaluation  for  the 
proposed  joint  estimation  approach  is  illustrated  in  Section  4. 
Finally,  conclusions  are  drawn  in  Section  5. 

2.  Modeling  the  lithium-ion  battery 

2.1.  The  dynamic  electrochemical-polarization  battery  model 

To  achieve  a  reliable  battery  state  estimation,  an  accurate  model 
must  be  built  first.  In  Refs.  [18,25],  seven  different  battery  models 
are  introduced  and  evaluated.  Due  to  the  remarkable  relaxation 
effect  of  the  lithium-ion  battery,  we  select  dynamic 
electrochemical-polarization  battery  model,  which  uses  an  elec¬ 
trochemical  battery  model  Nernst  model  to  replace  the  open  circuit 
voltage  (OCV)  part  of  the  Thevenin  model  making  SoC  available  as 
part  of  the  model  state  and  improving  the  model’s  accuracy.  The 
schematic  of  the  battery  model  is  shown  in  Fig.  1.  The  electrical 
behavior  of  the  proposed  model  can  be  expressed  by  Eq.  (1). 

1  &p  =  ~C^RpUp  +  Cp,L  (1) 

(  Ut  =  Hoc  -Up-  IlR0 

where  I<0,  /<j,  K2,  K3  and  I(4  are  five  constants  chosen  to  make  the 
model  well  fit  the  test  data,  R0  is  the  ohmic  resistance  and  Rp  is  the 
polarization  resistance.  The  polarization  capacitance  Cp  is  used  to 
describe  the  transient  dynamic  voltage  response  during  charging 
and  discharging.  Up  is  the  polarization  voltage  across  Cp,  Ut  is  the 
terminal  voltage.  Then  the  open  circuit  voltage  Uoc  can  be  described 
as  follows: 

Uoc  =  K0  +  K|  SoC  +  K2/S0C  +  K3ln  SoC  +  K4 ln(l  -  SoC)  (2) 


2.2.  The  parameters  identification  method 

For  battery  equivalent  circuit  models,  the  commonly  used 
parameters  identification  approach  is  multiple  linear  regression 
method  with  hybrid  pulse  power  characterization  data.  To  identify 
the  parameters  of  the  EP  model,  we  need  to  identify  the  Nemst 
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model’s  parameters  firstly.  The  parameters  of  Ki{i  =  0, 1, 2, 3, 4)  can 
be  calculated  by  linear  fitting  method  with  the  experiment  data  of 
the  OCV  vs.  the  SoC.  Then,  to  identify  the  other  parameters,  we 
firstly  have  a  discretization  process  for  dynamic  nonlinear  model 
and  then  we  should  build  a  regression  equation  for  the  discretiza¬ 
tion  system.  We  use  a  regression  equation  [18]  for  the  dynamic 
electrochemical-polarization  model  and  is  shown  in  Eq.  (3). 


2.3.2.  The  specific  hybrid  pulse  test 

The  battery  is  kept  in  the  thermal  chamber  and  the  tempera¬ 
ture  is  controlled  within  30  ±  5  °C.  Before  the  hybrid  pulse  test, 
the  available  capacity  test,  which  is  on  the  basis  of  the  standard  of 
Ref.  [26],  is  performed  to  measure  the  battery’s  maximum  avail¬ 
able  capacity.  The  maximum  available  capacity  Ca  of  the  battery  is 
the  number  of  ampere-hours  that  can  be  drawn  from  the  battery 


(3) 


where  Ipy  is  the  outflow  current  of  Cp  at  the  kth  sampling  time,  Ut,k 
is  the  terminal  voltage  at  the  kth  sampling  time,  1^  is  load  current 
at  the  kth  sampling  time  respectively  (assumed  positive  for 
discharge,  negative  for  charge).  The  At  is  the  sampling  interval.  The 
Uoc  is  determined  by  the  SoC— OCV  test.  And  the  time  constant  of 
polarization  (t  =  RPCP)  required  to  be  set  in  advance  for  regression 
operation,  and  for  different  time  constant,  different  discriminant 
coefficients  would  be  calculated;  then  an  optimum  time  constant 
value  will  be  achieved  through  finding  the  best  discriminant  coef¬ 
ficient  by  genetic  algorithm  [18]. 

2.3.  Battery  experiments 
2.3.1.  Test  bench 

The  schematic  of  the  built  test  bench  is  shown  in  Fig.  2,  which 
consists  of  an  Arbin  battery  test  system  BT2000  with  MITS  Pro  soft 
for  programming  the  BT2000,  a  thermal  chamber  for  environment 
control  and  a  host  computer.  The  BT2000  can  charge— discharge 
a  battery  according  to  the  designed  program  with  maximum 
voltage  of  60  V  and  maximum  charge/discharge  current  of  300  A 
with  three  scales  (5  A/50  A/300  A),  and  its  recorded  data  include 
current,  voltage,  temperature,  charge— discharge  Amp-hours  (Ah) 
and  Watt-hours  (Wh)  etc.  The  cell  voltage  also  can  be  measured  by 
the  auxiliary  channel,  and  its  measuring  range  is  0—5  V.  The 
measured  data  is  transmitted  to  the  host  computer  through  TCP/IP 
ports.  The  errors  of  the  Hall  current  and  voltage  sensors  are  less 
than  0.1%. 


Fig.  2.  The  schematic  of  the  battery  test  bench. 


at  the  standard  current  (C/3  rate),  starting  with  the  battery  is 
fully  charged,  which  maybe  different  with  its  nominal  value  due 
to  the  operation  environment.  For  the  fresh  LiMn204  battery, 
the  available  capacity  test  shows  that  its  actual  maximum  avail¬ 
able  capacity  is  35.4  Ah,  slightly  higher  than  35  Ah  of  the 
nominal  capacity.  Later,  an  OCV  test  is  carried  out  to  determine  the 
SoC-OCV  table. 

Then,  in  order  to  acquire  data  to  identify  the  model  parameters, 
a  specific  hybrid  pulse  test  procedure  is  conducted  at  0.1  SoC 
intervals  (constant  current  C/3  discharge  segments)  starting  from 
0.9  to  0.1  and  each  interval  followed  by  a  1-h  rest  to  allow  the 
battery  to  get  an  electrochemical  and  thermal  equilibrium  condi¬ 
tion  before  applying  the  next.  The  specific  hybrid  pulse  test  is 
similar  with  the  traditional  hybrid  pulse  power  characterization 
test,  but  the  specific  hybrid  pulse  test  uses  four  different  charge- 
discharge  currents  to  improve  the  applicability  of  the  EP  model 
under  a  more  broad  currents  operation  range  of  dynamic  driving 
cycles.  While  the  traditional  hybrid  pulse  power  characterization 
test  only  uses  1C  (discharge)  and  0.75C  (charge)  currents  to  build 
the  model,  and  the  model’s  parameters  error  will  be  arosed  due  to 
the  battery’s  current-dependent  relaxation  effect  and  coulombic 
efficiency  etc.  Considering  the  battery  current  operation  ranges  of 
the  battery  used  in  PHEVs  are  always  smaller  than  3C,  and  the 
maximum  current  is  less  than  4C,  we  choose  four  currents  (1C,  2C, 
3C  and  4C)  to  acquire  identification  data  set.  And  the  sampling 
interval  in  the  paper  is  1  s. 

The  specific  hybrid  pulse  test  results  are  shown  in  Fig.  3 
including  a  sampling  current  curve  in  Fig.  3(a),  a  sampling 
voltage  curve  under  SoC  =  0.8  in  Fig.  3(b),  the  current  profiles  in 
Fig.  3(c),  the  voltage  profiles  and  the  SoC  profiles  in  Fig.  3(d). 

2.3.3.  The  urban  dynamometer  driving  schedule  test 

The  urban  dynamometer  driving  schedule  (UDDS)  is  a  standard 
time  vs.  velocity  profile  for  urban  driving  vehicles  and  has  been 
widely  used  in  the  evaluation  of  model  accuracy  and  SoC  estimate 
of  battery  management  system  [8—10,19,27],  For  the  driving-cycle 
testing  of  PHEVs  batteries,  we  have  calculated  the  required  power 
from  the  powertrain  operating  with  the  UDDS  cycle.  Then  the 
required  power  is  reduced  as  a  certain  scale  and  then  we  use  the 
power  data  to  carry  out  the  battery  test.  For  the  UDDS  test, 
the  battery  SoC  is  controlled  in  the  range  of  0.15-0.9  according  to 
the  practical  requirements  of  the  electric  vehicles.  Due  to  the  hard 
determination  of  the  exact  SoC  value,  herein  we  determine  the 
initial  SoC  and  the  terminal  SoC  of  the  LiMn2C>4  lithium-ion  battery 
according  to  the  definition  of  SoC  with  a  standard  charging 
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Fig.  3.  The  specific  hybrid  pulse  test:  (a)  the  sampling  current  vs.  time  profile  of  one  cycle;  (b)  the  sampling  voltage  vs.  time  profile  of  one  cycle;  (c)  the  current  vs.  time  profiles;  (d) 
the  voltage  vs.  time  profiles  and  the  SoC  vs.  time  profiles. 


experiment  and  a  further  standard  discharging  experiment  after 
finishing  a  test,  so  the  initial  SoC  and  the  terminal  SoC  are  accurate. 
The  ampere-hour  counting  approach  is  used  to  calculate  the  SoC 
since  it  can  keep  track  of  the  accurate  SoC  with  an  accurate  initial 
SoC  and  a  compensation  of  the  coulombic  efficiency.  We  also 
improve  the  SoC  accuracy  with  a  revision  method  based  on  the 
accurate  terminal  SoC.  Considering  all  the  battery  experiments  are 
carried  out  in  a  temperature  chamber;  the  SoC  calculation  method 
is  feasible  with  an  acceptable  accuracy.  The  eight  consecutive  UDDS 
test  results  are  shown  in  Fig.  4  including  the  current  profiles, 
voltage  profiles  and  SoC  profiles. 


2.4.  The  parameters  identification  results  of  the  model 

The  Nemst  model’s  parameters  identification  results  with  the 
table  of  OCV  vs.  SoC  are  shown  in  Table  1  and  the  model’s  dynamic 
performance  parameters  results  are  shown  in  Table  2  with  the 
specific  hybrid  pulse  test. 

3.  SoC  and  peak  power  capability  joint  estimation  approach 

In  order  to  use  AEKF  for  the  SoC  and  peak  power  capability  joint 
estimation,  we  must  first  have  a  system  model  in  a  state-space  form 
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Fig.  4.  The  plots  of  the  UDDS  test:  (a)  the  sampling  current  vs.  time  profiles:  (b)  the  sampling  voltage  vs.  time  profiles;  (c)  the  sampling  SoC  vs.  time  profiles. 
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Table  1 

Identification  results  of  Nernst  model  parameters. 

K0  /Ci  K2  K3  K4 

3.8483  0.2441  -0.0124  0.0708  -0.0146 


[23],  Specifically,  we  assume  a  very  general  framework  for  discrete- 


time  lumped  dynamic  systems. 

Xfc+1  =  AX/{  +  Bufc  +  (,)k 

(4) 

Yr+i  =  CXfc+1  +  Du/,.  +  vk 

(5) 

We  know  that  an  AEKF-based  approach  will  be  used  to  estimate 
the  battery  SoC  given  noisy  measurements  [19,28,29],  To  do  so,  we 
require  a  discrete  state-space  model  similar  with  Eqs.  (4)  and  (5)  to 
describe  the  dynamics  of  the  battery  model.  Transform  Eq.  (1)  to 
a  discrete  system: 


UpM+1&  uPRexp(  +  4,fc+iKp(l  -  exp(  (8) 

Ut,k+i  uoc,fc+i  \,k+ iRo  -  l/p,W 

The  state  and  observation  function  of  the  discrete  system  is  as 
following: 


Y \  =  UtM,  uk  =  ILk,  A  =  ^6XP(  tO 

dUA  _  r_i  dtJoc(z)  | 


-exp(^)) 

v  At 


c  =  ^  .  -1 

9Xy=?:. 


dz 


(9) 


where  Xk  is  the  system  state  vector  at  the  /cth  sampling  time,  which 
represents  the  total  effect  of  system  inputs  uk  on  the  present  system 
operation,  such  as  SoC.  <ok  is  the  unmeasured  “process  noise”  that 
affects  the  system  state  and  vk  is  the  measurement  noise  which 
does  not  affect  the  system  state,  but  can  be  reflected  in  the  system 
output  estimates  Yk.  And  ajk  is  assumed  to  be  Gaussian  white  noise 
with  zero  mean  and  covariance  Q^;  is  assumed  to  be  Gaussian 
white  noise  with  zero  mean  and  covariance  R k.  The  matrices  A,  B,  C 
and  D  describe  the  dynamics  of  the  system,  and  are  time  varying 
and  determined  by  looking  up  for  the  parameters  table  vs.  SoC. 

3.1.  AEKF  based  SoC  estimation  approach 

To  carry  out  the  SoC  estimation  process,  According  to  the  SoC 
definition,  the  battery  SoC  is  the  ratio  of  the  remaining  capacity  to 
the  available  capacity: 

%  =  Zo  -  J  Vih,tdt/Ca  (6) 

o 

where  the  SoC  is  denoted  as  z,  zk  is  the  SoC  at  the  fcth  sampling  time, 
Zo  is  the  initial  SoC,  /L,t  is  the  instantaneous  load  current.  The  rji  is 
coulombic  efficiency,  which  is  the  function  of  the  current  and  the 
temperature  and  used  for  indicating  the  battery’s  capacity  loss  at 
different  operating  current.  To  describe  the  SoC  with  discrete  form, 
then  we  can  get  equation: 


where  d(J0C(z) /dz  =  K]  -  K2/z 2  +  K3/z  -  I<4/(  1  -  z)  calculated 
from  Eq.  (2). 

Then  we  can  achieve  AEKF-based  SoC  estimation  on  the  basis  of 
the  implementation  flowchart  of  the  AEKF  as  listed  in  Fig.  5. 

3.2.  The  EP  dynamic  model-based  peak  power  capability  estimation 
method 

3.2.1.  Instantaneous  peak  current  estimation  for  single  sampling 
interval 

To  build  an  EP  model-based  dynamic  multi-parameter  method 
for  peak  power  capability  estimation  of  lithium-ion  battery,  we 
must  first  have  a  discrete  battery  model  written  as  follows  on  the 
basis  of  Eq.  (8). 

Ot,fc+i  =  Hoc (zfc+1 )  -  Up  kexp 

-Wi(Ko  +  Hp(l-exp(^))) 

However,  the  peak  currents  cannot  be  solved  directly  from  the 
maximum  current  kik  since  zk+ i  itself  is  a  function  of  the  current 
Iiy  and  OCV  is  a  nonlinear  function  of  SoC  [2],  To  avoid  this 
problem,  the  Taylor-series  expansion  is  employed  to  linearize  the 
equation  and  to  solve  the  approximated  values  for  the  peak 
currents.  The  Taylor-series  expansion  equation  for  OCV  is  as 
follows: 


Zk  =  zk  i  -  vA^t/c, 


(7) 


Parameters  identification  results  with  SoC  =  0.6,  0.7,  0.S 


SoC  Cp/F  Rp/mO  Rp/mQ  t/s 

"06  29,029  1.142  1.914  33.16 

0.7  26,484  1.335  1.925  35.36 

0.8  29,250  1.045  1.936  30.57 


Uoc(zk+1)  =  Uoc(zk-ILM1r^j 


+R1  (zkA,k+ 1"^“) 

(ID 


Considering  the  first-order  residual  i?i(  )  being  too  small  to 
affect  the  OCV  at  the  next  sampling  time  At,  since  the  SoC’s 
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Note:  Hk  is  the  innovation  covariance  matrix  based  on  the  innovation  sequence  inside  a  moving  estimation 
window  of  size  M.  Qk  and  Rk  for  the  Q  and  R  at  the  Ath  sampling  time  respectively.Where  Kkis  the  kalman 
gain  matrix;  ek  is  defined  as  the  difference  between  the  measurement  and  the  observation  CXk  +  Duk ,  , 

Xk  and  Xk  are  for  the  priori  estimate  before  the  measurement  is  taken  into  account  and  the  posteriori 
estimate  after  the  measurement  is  taken  into  account  respectively. 4.1  is  the  initial  time  of  calculation. 


Fig.  5.  An  implementation  flowchart  of  the  AEKF  algorithm. 


variation  per  sampling  time  is  very  small  [19,27],  the  Ri(  )  can  be  where  /^EP  and  /miiP  are  the  minimum  charge  currents  and 
viewed  as  Ri  ( ■ )  =  0.  Then:  maximum  discharge  currents  under  the  EP  dynamic  model-based 

method. 


V  a  I  lz=zk  (12) 

+  Rp  (1  -  exp  )  )  -  Up,kexp 

Therefore,  we  can  get  the  peak  current  estimates  when  takes  the 
considerations  of  terminal  voltage  design  limits  —  Lit, min  and  Ut.max, 
and  the  computational  equation  is  as  follows: 


Uoc(zk)  -  UpMexp(-^j  -  Utmin 

Uoc(zk)  -  Upkexp(-^)  -  Ut,max 
,chg,EP  _  _  V  t  / _ 

-  -~|a  +  Bp(,-eXp(4))+B. 


3.2.2.  Continuance  peak  current  estimation  for  multi  sampling 
intervals 

The  peak  current  estimation  approach  calculated  by  Eq.  (13) 
shows  that  the  estimates  have  a  big  difference  with  different 
sampling  interval  At;  and  for  electric  vehicles  application,  we  usually 
require  to  restrict  the  peak  currents  for  some  special  conditions, 
such  as  acceleration,  climbing,  braking  etc.,  for  this  respect, 
a  continuance  peak  power  estimation  approach  for  multi  sampling 
intervals,  which  is  on  the  basis  of  the  instantaneous  peak  current 
estimation  approach  for  single  sampling  interval,  is  proposed  here. 

The  input  of  the  system  is  assumed  constant  between  the  kth 
sampling  time  tk  and  the  (k  +  L)th  sampling  time  tk+L,  which  is 
expressed  by:  uk+k  =  uk  with  the  sampling  intervals  is  L  x  At.  Then, 
we  can  use  the  EP  model  to  predict  voltage  at  the  (fe  +  L)th  sampling 
time  into  the  future  by: 

f  Xfc  I  L  =  AX/c+i-l  +  Bu/df.  1 

l  Yk+L  =  CXk+L  +  Du/;  /  1 


(14) 
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Since  Eq.  (1)  should  be  written  in  a  linear  discretization  form  as 
Eq.  (14),  therefore,  for  input  constant  over  the  entire  prediction 
horizon  (from  time  tk  to  tk+L),  we  have: 


To  find  the  minimum  charging  current  the  output  of  the 

system  should  satisfy  the  following  equation: 


Uoc(%n) - Up,k (exp (~^j ^  -Cin,f  +RP^l-exp ^ 


To  find  the  maximum  discharging  current  ijj^p,  the  output  of 
the  system  should  satisfy  the  following  equation: 

Uoc(Zfc+t)  -  Up  k  (exp -;maxEL  ^o+Rp  (l -exp(-^)) 

O’) 

Then,  we  can  get  the  computational  approach  to  predict  peak  current 
at  the  kth  sampling  time  into  the  ( k  +  L)th  sampling  time  by: 


battery  will  be  over-discharged;  on  the  other  hand,  we  should 
control  its  charge  currents  and  maximize  its  discharge  currents 
when  the  actual  SoC  is  close  to  its  maximum  design  limits  zmax, 
otherwise,  the  battery  will  be  over-charged.  This  is  the  key  prin¬ 
ciple  for  the  peak  current  estimation  approach  with  the  SoC-limited 
method,  and  computational  equation  and  relevant  description  is 
shown  as  follows,  which  is  derived  from  Eq.  (7)  when  takes  the 
consideration  of  SoC  design  limits. 


jChg,soc 

Z-k  ~  ^max 

“  %EAt/Ca 

/^°c 

_Zk~  Zmin 

vMt/cs 

(19) 


where  and  are  the  minimum  charge  current  and 

maximum  discharge  current  between  the  L  x  At  sampling  intervals 
under  the  SoC  design  limits.  Once  the  current  design  limit  is  calcu¬ 
lated,  the  peak  currents  with  all  limits  enforced  are  calculated  as: 


f«.- 

\ Cn  -  . Cm"  Cm) 

where  /max  and  fmjn  are  design  limits  for  the  maximum  discharge 
current  and  minimum  charge  current  respectively.  and  /[jj|x  are 
the  minimum  charge  current  and  maximum  discharge  current 
respectively  with  all  limits  enforced.  Then  the  peak  power  capa¬ 
bility  estimation  based  on  the  EP  dynamic  model-based  can  be 
calculated  as  follows: 


ydis.EP 


jThg.EP 


Uoc(zk)  -  tlP:j<(exp(  f))Ut.min 

Uoc(Zfc)-Up,k(exp(^))  -Otmax 

^d^|Zt  +  Rp(l_eXp(^))|(exp(^))* i-1 * *-V,0 


(18) 


When  the  L  is  set  as  one,  Eq.  (18)  is  the  same  as  Eq.  (13). 

3.2.3.  Peak  power  capability  estimation  method 

To  make  the  battery’s  performance  more  safe  and  reasonable 
when  the  actual  SoC  is  close  to  its  design  limits,  we  should  control 
its  discharge  currents  and  maximize  its  charge  currents  when  the 
actual  SoC  is  close  to  its  minimum  design  limits  zm jn,  otherwise,  the 


/  PZ  =  max(pmin, UtkaCn) 

l  Pfjj!x  =  min(PmaX;[/t,k+^x) 

Then  we  can  get  the  computational  approach  for  peak  power 

estimation  when  we  take  the  terminal  voltage  calculation  equation 

shown  in  Eq.  (10). 


pSSn~max|Pmin,  ^oc(zfc+i)  -  Up,fc(exp(  (r0  +  Rp (l  -  exp ^  (exp (^) y  1 

P^x  «  min  |pmax,  (u0C(zk+L)  -  Up.fc  (exp  (=^)  )  ‘-ig**  +  RP  (l  -  exp  (^)  )  E*  (exp  (~f)  pJJ 


(22) 
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where  Pmax,  Pmin  are  the  battery’s  power  design  limits,  Pmax  denotes 
the  peak  discharge  power  and  Pmjn  denotes  the  peak  charge  power 
of  its  design  limits  respectively. 

3.3.  State  of  charge  and  peak  power  capability  joint  estimation 
algorithm 

In  order  to  use  AEKF-based  methods  for  the  SoC  and  peak  power 
capability  joint  estimation,  we  must  first  have  a  complete  sche¬ 
matic  of  the  online  estimation  and  computational  process.  There¬ 
fore,  a  very  general  framework  for  online  joint  estimation,  which 
named  six  steps  joint  estimation  method,  is  shown  in  Fig.  6.The 
relevant  processes  are  summarized  as  follows. 

Step  l :  Initialization  —  when  given  a  random  initial  X  matrix,  we 
can  get  the  model’s  parameters  through  looking  up  for  the 
parameters  tables  and  which  is  ready  for  AEKF-based  online 
estimation  approach.  This  initialization  will  be  finished  after  the 
program  code  is  loaded  on  the  equipment. 

Step  2:  Start  —  when  an  input  signals  of  the  system,  load  current 
or  required  power  is  given,  the  charge-discharge  current  is 
loaded  on  the  LiMn2C>4  lithium-ion  battery  and  the  battery 
model  simultaneously. 

Step  3:  Correction  and  SoC  estimation  —  The  terminal  voltage 
error  between  the  observer  and  the  experimental  data  is 
adaptively  reduced  by  updating  the  gain  of  the  AEKF.  And  next 
the  updated  gain  is  used  to  compensate  for  the  state  estimation 
error.  The  SoC  estimate  is  then  fed  back  to  update  the  parame¬ 
ters  of  the  battery  model  for  the  SoC  estimation  at  the  next 
sampling  time.  And  an  accurate  estimate  for  SoC  and  Up  will  be 
reached  quickly. 

Step  4:  Dynamic  peak  current  estimation  -  Then,  we  require  to 
calculate  the  peak  currents  with  EP  dynamic  method  and 
designed  SoC-limited  method,  which  will  be  performed  by  Eqs. 
(18)  and  (19). 

Step  5:  Dynamic  multi-parameter  peak  power  capability  estima¬ 
tion  method  -  When  we  take  the  consideration  of  the  design 
limits  of  current,  voltage  and  power,  we  can  achieve  the  peak 
power  estimates  with  Eqs.  (20)  and  (22). 

Step  6:  End  -  Lastly,  we  can  achieve  the  joint  estimation  for  the 
SoC  and  peak  power  capability. 


Table  3 

Parameters  limits  for  peak  power  capability  estimation. 

Parameters  Maximum  value  Minimum  value 

SoC  (Zmax,  Zmin)  09  015 

Ut(Utnax,Utmin)/V  4.25  3.0 

3l  (3max,  3min)/A  350  -175 

P  (Pmax.  Pmin)/W  1500  -750 


4.  Verification  and  discussion 

This  section  will  verily  and  evaluate  the  AEKF-based  SoC  and 
peak  power  capability  joint  estimation  approach  with  different 
initial  SoC  under  the  UDDS  test  data.  First,  we  should  define  the  SoC 
design  limits  of  the  PHEVs  batteries  except  for  such  design  limits  as 
voltage,  power  and  current  defined  by  the  battery  manufacturers.  It 
is  obvious  that  for  different  control  strategies,  the  SoC  limits  maybe 
different.  An  example  of  parameter  limits  for  peak  power  capability 
estimation  is  listed  in  Table  3. 

Where  the  peak  current  and  peak  power  design  limits  are 
restricted  for  continuous  power,  therefore  for  a  short  time  which 
less  than  60  s,  the  peak  currents  and  powers  will  be  much  bigger 
than  the  values  present  in  Table  3  and  maybe  700/— 350  A  and 
3000/-1500  W. 

4.1.  Verification  for  joint  estimation  method 

The  AEKF-based  SoC  estimation  approach  was  conducted  and 
the  results  are  shown  in  Fig.  7  with  the  initial  parameters  as 
following: 

x0  =  [4.1  0.9 ]T,  P0  =  [J  °],  Q°  =  [J  R0  =  1  (23) 

From  Fig.  7(a),  we  can  see  the  estimated  terminal  voltage  tracks 
the  experiment  profiles  well,  and  the  details  are  shown  in  Fig.  7(b). 
From  Fig.  7(b),  it  indicates  that  the  terminal  voltage  error  is  less 
than  2%  of  its  nominal  voltage,  which  is  because  the  AEKF-based 
algorithm  can  precisely  estimate  the  voltage  and  timely  adjust 
the  Kalman  gain  matrix  I</,  according  to  the  terminal  voltage  error 
between  the  measured  data  and  the  estimated  values.  From 


SoC  estimation  Peak  power  capability  estimation 


Fig.  6.  The  flowchart  of  the  joint  estimation  under  the  AEKF-based  adaptive  observer. 
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Fig.  7.  AEKF-based  SoC  estimation  results  with  an  accurate  initial  SoC:  (a)  terminal  voltage;  (b)  voltage  estimation  error;  (c)  SoC;  (d)  SoC  estimation  error. 


Fig.  7(c),  we  can  see  the  estimated  SoCs  track  the  true  SoCs  well, 
Fig.  7(d)  shows  that  the  range  of  the  SoC  estimation  error  is  nearly 
between  -0.01  and  +0.01,  which  is  quite  accurate  for  SoC 
estimation. 

And  we  also  have  calculated  the  real-time  peak  power  capability 
with  the  six  steps  methods  for  15  s  30  s  and  60  s,  the  results  are 
shown  in  Fig.  8. 

From  Fig.  8,  we  can  see  the  more  sampling  intervals  are 
required,  the  smaller  of  the  absolute  values  of  its  peak  currents  and 
peak  powers  will  be.  Therefore,  the  instantaneous  peak  powers 
estimation  are  not  suitable  for  longtime  application  or  else  the 
battery  will  be  over-charged  or  over-discharged  and  damaged. 


Additionally,  the  EP  model-based  method  can  simulate  the 
dynamic  performance  of  the  battery  accurately  and  the  dynamic 
power  capability  estimation  results  are  well  for  practical  use. 
However,  the  above  estimation  for  SoC  and  peak  power  capability  is 
based  on  an  accurate  initial  SoC,  and  an  accurate  SoC  estimation 
depends  on  two  aspects  according  to  the  definition  of  SoC  given  by 
Eq.  (6),  one  is  the  calculation  of  accumulated  capacity  consumption, 
another  and  most  important  is  the  robust  performance  for  initial 
error  of  the  SoC.  The  above  results  and  analysis  show  that  the  SoC 
estimation  error  under  the  precise  initial  SoC  is  less  than  0.02, 
which  is  quite  desirable.  In  order  to  determine  whether  the  joint 
estimation  with  AEKF  approach  can  effectively  solve  the  initial 


Fig.  8.  AEKF-b, 
charge  power. 


peak  power  capability  estimatic 


i  initial  SoCs:  (a)  peak  discharge  i 


t;  (b)  peak  charge  current;  (c)  peak  discharge  power;  (d)  peak 
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Fig.  9.  Robust  performance  evaluation  on  SoC  estimation  results:  (a)  SoC  estimation  with  5o(n  =  0.95  and  SoC0  =  0.50;  (b)  zoom  figure  for  (a);  (c)  SoC  estimation  error  with 
SoC0  =  0.95  and  SoQ,  =  0.50;  (d)  zoom  figure  for  (c). 


error  of  the  SoC,  we  carry  out  a  further  analysis  for  evaluating  its 
robust  performance. 

4.2.  Evaluation  on  the  robustness  of  the  joint  estimation 
performance 

Two  different  SoC  initial  values,  0.95  and  0.50,  are  preset  and  the 
corresponding  SoC  estimations  are  performed  based  on  the  UDDS 
cycles.  The  robust  performance  evaluation  results  for  SoC  are 
shown  in  Fig.  9. 

From  Fig.  9(a),  we  can  see  the  SoC  estimation  results  with  the 
two  different  initial  SoC  are  all  almost  the  same  as  the  true  data,  the 
zoom  figure  strengthens  the  description  of  the  robustness  perfor¬ 
mance  of  the  AEKF-based  SoC  estimation  approach.  Fig.  9(b)  indi¬ 
cates  the  proposed  method  ensure  the  robustness  of  SoC  estimation 
after  ten  sampling  intervals.  To  make  the  evaluation  more  clarity, 
we  plot  the  SoC  estimation  error  in  Fig.  9(c)  and  (d),  which  indicates 
that  the  AEKF-based  SoC  estimation  approach  precisely  estimates 
the  SoC  and  timely  modulates  to  adjust  the  Kalman  gain  matrix 
according  to  the  terminal  voltage  error.  And  the  SoC  estimation 
error  would  be  reduced  to  0.02  after  10  sampling  intervals,  which  is 
well  suitable  for  PHEVs  battery  management  system  application. 
The  precision  of  the  SoC  estimation  ensures  the  reliable  and 
robustness  performance  of  the  peak  power  capability  estimation 


results.  We  use  the  peak  discharge  current  estimates  as  an  example 
to  evaluate  the  robustness  performance  of  the  peak  power  esti¬ 
mation  capability.  Fig.  10  is  the  peak  discharge  currents  estimation 
with  different  initial  SoC  values  and  the  calculation  sampling 
intervals  is  30  s,  the  comparison  profiles  of  the  peak  discharge 
current  estimation  results  with  three  initial  SoC  values  and  the  true 
SoC  are  plotted  in  Fig.  10(a),  and  zoom  figure  is  in  Fig.  10(b). 

Due  to  the  peak  power  capability  estimation  effect  has  a  direct 
correction  with  the  SoC  estimation  result  from  the  calculation 
approach  shown  in  Fig.  6,  the  accurate  SoC  estimation  has  been 
a  strong  guarantee  for  the  robustness  performance  of  the  peak 
power  capability  estimation.  Fig.  10(a)  shows  that  the  peak  current 
estimates  are  almost  the  same  for  three  SoC  estimation  results  and 
the  true  SoC  data,  and  the  zoom  figure  shows  that  the  estimation 
results  can  converge  to  the  stable  value  within  10  sampling  inter¬ 
vals  especially  with  error  initial  SoC. 

Therefore,  the  proposed  AEKF-based  joint  estimation  approach 
for  SoC  and  peak  power  capability  can  efficiently  estimate  the  SoC 
and  peak  power  capability  together  accurately  especially  with  big 
initial  SoC  error,  and  the  maximum  SoC  estimation  error  is  less  than 
2%  after  10  sampling  intervals.  In  addition,  it  is  obvious  that  accu¬ 
rate  SoC  and  available  peak  power  capability  estimates  will  benefit 
the  Vehicle  to  Grid  (V2G)  application,  and  improve  the  control 
efficiency  of  the  peak  load  regulation  of  the  grid. 


Time/mm  Time/sec 


Fig.  10.  Robust  performance  evaluation  on  peak  power  estimation  results:  (a)  peak  discharge  current  estimates  of  three  initial  SoC  values  and 


e;  (b)  zoom  figure  for  (a). 
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5.  Conclusions 

Based  on  the  above  analysis,  the  main  concluding  remarks  can 
be  made  below: 

(1)  To  model  the  dynamic  performance  of  Lithium-ion  battery 
accurately,  the  dynamic  electrochemical  polarization  battery 
model  is  employed  for  state  estimations:  the  model  uses  the 
Nernst  model  to  build  a  relationship  between  the  SoC  and  the 
open  circuit  voltage  characteristic  of  the  battery,  therefore 
the  model  not  only  has  the  advantages  of  simulation  accuracy 
for  the  terminal  voltage,  but  also  having  a  contribution  to 
improve  state  estimation  accuracy. 

(2)  To  get  accurate  model’s  parameters  and  improve  the  dynamic 
performance  simulation  precision  of  battery  model,  we  rede¬ 
fine  the  hybrid  power  pulse  characteristic  test  with  four 
different  charge-discharge  currents;  and  then  we  build 
a  complete  parameters  identification  process  for  the  EP  model. 

(3)  To  achieve  a  robust  SoC  estimation,  the  AEKF  approach  is 
employed  on  the  basis  of  the  electrochemical  polarization 
battery  model. 

(4)  To  achieve  a  continuous  peak  power  requirement  of  the  PHEVs, 
the  continuous  peak  power  capability  estimation  approach  for 
multi-sampling  intervals  is  proposed,  and  the  corresponding 
peak  current  and  peak  power  calculation  methods  are  also 
built. 

(5)  An  AEKF-based  six-step  joint  estimation  approach  for  SoC  and 
peak  power  capability  is  proposed  and  its  general  operation 
flowchart  is  built.  The  simulation  and  experiment  results 
indicate  that  the  proposed  approach  not  only  has  the  advan¬ 
tages  of  online  estimating  the  SoC  accurately  and  robustly,  and 
its  peak  error  is  less  than  2%  especially  with  big  initial  SoC  error 
after  10  sampling  intervals,  but  also  can  predict  the  peak  power 
capability  reliably  and  robustly. 

The  future  work  will  focus  on  the  joint  estimation  approach 
with  the  battery  pack  and  online  parameters  identification  method 
for  the  dynamic  battery  model,  and  the  systematic  validation  test 
scheme  for  available  peak  power  capability  estimation. 
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